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Abstract 

A class of three-dimensional initial data characterized by uniformly 
large vorticity is considered for the 3D incompressible Euler equations in 
bounded cylindrical domains. The fast singular oscillating limits of the 
3D Euler equations are investigated for parametrically resonant cylin- 
ders. Resonances of fast oscillating swirling Beltrami waves deplete the 
Euler nonlinearity. These waves are exact solutions of the 3D Euler 
equations. We construct the 3D resonant Euler systems; the latter are 
countable uncoupled and coupled 5*0(3; C) and SO(3; R) rigid body 
systems. They conserve both energy and helicity. The 3D resonant Eu- 
ler systems are vested with bursting dynamics, where the ratio of the 
enstrophy at time t = t* to the enstrophy at t = of some remarkable 
orbits becomes very large for very small times t*; similarly for higher 
norms H s , s > 2. These orbits are topologically close to homoclinic 
cycles. For the time intervals where H s norms, s > 7/2 of the limit 
resonant orbits do not blow up, we prove that the full 3D Euler equa- 
tions possess smooth solutions close to the resonant orbits uniformly in 
strong norms. 
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1 Introduction 



The issues of blowup of smooth solutions and finite time singularities of 
the vorticity field for 3D incompressible Euler equations are still a major 
open problem. The Cauchy problem in 3D bounded axisymmetric cylindrical 
domains is attracting considerable attention: with bounded, smooth, non- 
axisymmetric 3D initial data, under the constraints of conservation of bounded 
energy, can the vorticity field blow up in finite time? Outstanding numeri- 
cal claims for this have recently been disproven [Kej . [Houlj . |Hou2] . The 
classical analytical criterion of Beale-Kato-Majda |B-K-M| for non-blow up 
in finite time requires the time integrability of the L°° norm of the vorticity. 
DiPcrna and Lions [Oj have given examples of global weak solutions of the 
3D Euler equations which are smooth (hence unique) if the initial conditions 
are smooth (specifically in W 1,P (_D), p > 1). However, these flows are really 
2-Dimensional in x±,X2, 3-components flows, independent from the third co- 
ordinate x 3 . Their examples |DiPe-Lij show that solutions (even smooth ones) 
of the 3D Euler equations cannot be estimated in W 1,p for 1 < p < oo on any 
time interval (0, T) if the initial data are only assumed to be bounded in W 1,p . 
Classical local existence theorems in 3D bounded or periodic domains by Kato 
|Kaj . Bourguignon-Brezis [Bou-Br] and Yudovich [Yulj . |Yu2| require some 
minimal smoothness for the initial conditions (IC), e.g., in H S (D), s > |. 
The classical formulation for the Euler equations is 

d t V + (V ■ V)V = -Vp, V • V = 0, (1.1) 
V-N = 0on<9L>, (1.2) 

where dD is the boundary of a bounded, connected domain D, N the normal 
to dD, V(t,y) = (Vi,V 2 ,V 3 ) the velocity field, y = (yi,y 2 ,2/3), and p is the 
pressure. 

The equivalent Lame form [Ar-Khc 

d t V + curlV x V + V (p+ i|V| 2 ^ = 0, (1.3) 
V-V = 0, (1.4) 

d t u + curl{u X V) = 0, (1.5a) 
uj = curlV, (1.5b) 

implies conservation of Energy: 



D 



2 

The helicity Hel(t) |Ar-Khej . [Mofj . is conserved 



E(t) = ^ I \V{t,y)\ 2 dy. (1.6) 



Hel{t) = / V -u>dy, (1.7) 
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for D = R 3 and when D is a periodic lattice. Helicity is also conserved for 
cylindrical domains, provided that u>-N — on the cylinder's lateral boundary 
at t = (see |M-N-B-Gp . 

From the theoretical point of view, the principal difficulty in the analysis 
of 3D Euler equations is due to the presence of the vortex stretching term 
(u> ■ V)V in the vorticity equation (|1 .5a[) . The equations (| 1 . 3[) and (|1.5ap are 
equivalent to: 

d t u + [u;,V] =0, (1.8) 

where [a, b] = curl (a x b) is the commutator in the infinite dimensional Lie 
algebra of divergence-free vector fields |Ar-Khej . This point of view has led to 
celebrated developments in Topological Methods in Hydrodynamics |Ar-Khe] , 
|Mofj . The striking analogy between the Euler equations for hydrodynamics 
and the Euler equations for a rigid body (the latter associated to the Lie 
Algebra of the Lie group 5*0(3, R)) had already been pointed out by Moreau 
|Morlj : Moreau was the first to demonstrate conservation of Helicity (1961) 
|Mor2]. This has led to extensive speculations to what extent/in what cases 
are the solutions of the 3D Euler equations "close" to those of coupled 3D rigid 
body equations in some asymptotic sense. Recall that the Euler equations for 
a rigid body in R 3 is: 

m ( + a; x m = 0, m - Alj, (1.9a) 
m t + [w,m] = 0, (1.9b) 

where m is the vector of angular momentum relative to the body, ui the 
angular velocity in the body and A the inertia operator [Arlj , |Ar-Khej . 

The Russian school of Gledzer, Dolzhansky, Obukhov |G-D-Oj and Vishik 
[Vish] has extensively investigated dynamical systems of hydrodynamic type 
and their applications. They have considered hydrodynamical models built 
upon generalized rigid body systems in SO(n, R), following Manakhov [Man] . 
Inspired by turbulence physics, they have investigated "shell" dynamical sys- 
tems modeling turbulence cascades; albeit such systems are flawed as they 
only preserve energy, not helicity. To address this, they have constructed and 
studied in depth n-dimensional dynamical systems with quadratic homoge- 
neous nonlinearities and two quadratic first integrals F\ , F2 ■ Such systems 
can be written using sums of Poisson brackets: 

dxh = 1 v (dF L m_dF L m 

dt 2 ^ e Pii - t,l \dx^dx^ cV3 dx^ 

where constants Pi A ...i n are antisymmetric in i^, i n . 

A simple version of such a quadratic hydrodynamic system was introduced 
by Gledzer [Gil] in 1973. A deep open issue of the work by the Gledzer- 
Obukhov school is whether there exist indeed classes of I.C. for the 3D Cauchy 
Euler problem (|1.1[) for which solutions are actually asymptotically close in 
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strong norm, on arbitrary large time intervals to solutions of such hydro- 
dynamic systems, with conservation of both energy and helicity. Another 
unresolved issue is the blowup or global regularity for the "enstrophy" of such 
systems when their dimension n — > oo. 

This article reviews some current new results of a research program in 
the spirit of the Gledzer-Obukhov school; this program builds-up on the re- 
sults of [M-N-B-G for 3D Euler in bounded cylindrical domains. Following 
the original approach of |B-M-Nlj - [B-M-N4j in periodic domains, |M-N-B-G| 
prove the non blowup of the 3D incompressible Euler equations for a class 
of three-dimensional initial data characterized by uniformly large vorticity in 
bounded cylindrical domains. There arc no conditional assumptions on the 
properties of solutions at later times, nor are the global solutions close to 
some 2D manifold. The initial vortex stretching is large. The approach of 
proving regularity is based on investigation of fast singular oscillating limits 
and nonlinear averaging methods in the context of almost periodic functions 
|Bo-Mij . [Besj . |Cor| . Harmonic analysis tools based on curl eigenfunctions 
and eigenvalues are crucial. One establishes the global regularity of the 3D 
limit resonant Euler equations without any restriction on the size of 3D initial 
data. The resonant Euler equations are characterized by a depleted nonlin- 
earity. After establishing strong convergence to the limit resonant equations, 
one bootstraps this into the regularity on arbitrary large time intervals of the 
solutions of 3D Euler Equations with weakly aligned uniformly large vorticity 
at t = 0. M- N-B-G] theorems hold for generic cylindrical domains, for a set 
of height/radius ratios of full Lebesgue measure. For such cylinders, the 3D 
limit resonant Euler equations are restricted to two-wave resonances of the 
vorticity waves and are vested with an infinite countable number of new con- 
servation laws. The latter are adiabatic invariants for the original 3D Euler 
equations. 

Three- wave resonances exist for a nonempty countable set of h/R (h 
height, R radius of the cylinder) and moreover accumulate in the limit of 
vanishingly small vertical (axial) scales. This is akin to Arnold tongues [Ar2j 
for the Mathieu-Hill equations and raises nontrivial issues of possible sin- 
gularities/lack thereof for dynamics ruled by infinitely many resonant triads 
at vanishingly small axial scales. In such a context, the 3D resonant Euler 
equations do conserve the energy and helicity of the field. 

In this review, we consider cylindrical domains with parametric resonances 
in h/R and investigate in depth the structure and dynamics of 3D resonant 
Euler systems. These parametric resonances in h/R are proven to be non- 
empty. Solutions to Euler equations with uniformly large initial vorticity are 
expanded along a full complete basis of elementary swirling waves (T 2 in time). 
Each such quasiperiodic, dispersive vorticity wave is a quasiperiodic Beltrami 
flow; these are exact solutions of 3D Euler equations with vorticity parallel 
to velocity. There are no Galerkin-like truncations in the decomposition of 
the full 3D Euler field. The Euler equations, restricted to resonant triplets of 
these dispersive Beltrami waves, determine the "resonant Euler systems" . The 
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basic "building block" of these (a priori oo-dimensional) systems are proven 
to be SO(3; C) and .SO^R) rigid body systems: 

Uk + (A m — \ n )U m U n = 

U m + (A„ - X k )U n U k = (1.11) 

Un + (Afe - \ m )U k U m = 

These A's are eigenvalues of the curl operator in the cylinder, curl$>^ = 
±A„$„; the curl eigenfunctions are steady elementary Beltrami flows, and 
the dispersive Beltrami waves oscillate with the frequencies ±2§j x 1 ' n 3 ver ~ 
tical wave number (vertical shear), < e < 1. Physicists |Ch-Ch-Ey-H have 



computationally demonstrated the physical impact of the polarization of Bel- 
trami modes $ ± on intermittency in the joint cascade of energy and helicity 
in turbulence. 

Another "building block" for resonant Euler systems is a pair of SO(3; C) 
or 50(3; R) rigid bodies coupled via a common principal axis of inertia/mo- 
ment of inertia: 

(1.12a) 
(1.12b) 

A m )ra^a m (1.12c) 
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(1.12d) 
(1.12e) 

where T and T are parameters in R defined in Theorem 14.101 Both reso- 
nant systems (jl.lip and (|1.12|) conserve energy and helicity. We prove that 
the dynamics of these resonant systems admit equivariant families of homo- 
clinic cycles connecting hyperbolic critical points. We demonstrate bursting 
dynamics: the ratio 

||u(t)||W||u(0)||i, s , *>1 

can burst arbitrarily large on arbitrarily small times, for properly chosen para- 
metric domain resonances h/R. Here 

IKt)|| 2 H a=5> n ) 2 > n (t)| 2 . (1-13) 
a 

The case s = 1 is the enstrophy. The "bursting" orbits are topologically close 
to the homoclinic cycles. 

Are such dynamics for the resonant systems relevant to the full 3D Euler 
equations (|l.l| - (|1.8p ? The answer lies in the following crucial "shadowing" 
Theorem 12.101 Given the same initial conditions, given the maximal time 
interval < t < T m where the resonant orbits of the resonant Euler equations 
do not blow up, then the strong norm H s of the difference between the exact 
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Eulcr orbit and the resonant orbit is uniformly small on < t < T m , provided 
that the vorticity of the I.C. is large enough. Paradoxically, the larger the 
vortex streching of the I.C, the better the uniform approximation. This deep 
result is based on cancellation of fast oscillations in strong norms, in the 
context of almost periodic functions of time with values in Banach spaces 
(Section 4 of [M-N-B-G ). It includes uniform approximation in the spaces 
H s , s > 5/2. For instance, given a quasiperiodic orbit on some time torus T' 
for the resonant Euler systems, the exact solutions to the Euler equations will 
remain e-close to the resonant quasiperiodic orbit on a time interval < t < 
maxTj, 1 < i < I, Tj elementary periods, for large enough initial vorticity. If 
orbits of the resonant Euler systems admit bursting dynamics in the strong 
norms H s , s > 7/2, so do some exact solutions of the full 3D Euler equations, 
for properly chosen parametrically resonant cylinders. 

2 Vorticity waves and resonances of elemen- 
tary swirling flows 

We study initial value problem for the three-dimensional Euler equations with 
initial data characterized by uniformly large vorticity: 

d t V + (V • V)V = -Vp, V • V = 0, (2.1) 
V(t, y)\ t=0 = V(0) = V (y) + ^e 3 xy (2.2) 

where y = (yi, y 2 , 1/3), V(i, y) = (Vi, V%, V3) is the velocity field and p is the 
pressure. In Eqs. e 3 denotes the vertical unit vector and CI is a constant 
parameter. The field Vo(y) depends on three variables y\, 7/2 and z/3. Since 
curl(-^e 3 X y) = fle^, the vorticity vector at initial time t = is 

curlV(0, y) = curlV {y) + fle 3 , (2.3) 

and the initial vorticity has a large component weakly aligned along e3, when 
f2 >> 1. These are fully three-dimensional large initial data with large initial 
3D vortex stretching. We denote by the usual Sobolev space of solenoidal 
vector fields. 
The base flow 

V s (y) = ^e 3 x y, curlV s (y) = fle 3 (2.4) 

is called a steady swirling flow and is a steady state solution (|l.ip - (|1.4[) . as 
curl{VLe 3 x V s (y)) = 0. In (|2.2p and (|2.3p . we consider I.C. which are an arbi- 
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trary (not small) perturbation of the base swirling flow ~V s (y) and introduce 

V(t,y) = ^e 3 xy + V(t,y), (2.5) 

curlV(t, y) = fle 3 + curlV{t, y), (2.6) 
d t V + curlV x V + Qe 3 x V + curlV x V s (y) + Vp' = 0, V • V = 0, (2.7) 
V(t, y )| t =o=Vo(»). (2.8) 

Eqs. (|2.ip and (|2.7[) are studied in cylindrical domains 

C = {(t/x, 2/2,2/3) G R 3 : < y 3 < 2n/ a , y\ + y 2 < R 2 } (2.9) 

where a and i? are positive real numbers. If h is the height of the cylinder, 
a = 2ir/h. Let 

r = {(2/i,y 2 ,2/3)eR 3 :0<2/ 3 <27r/a, y 2 + y 2 = R 2 }. (2.10) 

Without loss of generality, we can assume that R = 1. Eqs. (|2.1[) are consid- 
ered with periodic boundary conditions in y 3 

V(yi,2/2,y3) = V(yi, ?/2,y3 + 27r/a) (2.11) 

and vanishing normal component of velocity on T 

V ■ N = V ■ N = on T; (2.12) 

where N is the normal vector to T. From the invariance of 3D Eulcr equations 
under the symmetry y 3 — * —3/3, V\ — > Vi, V2 — > V2, V3 — > — V3, all results in 
this article extend to cylindrical domains bounded by two horizontal plates. 
Then the boundary conditions in the vertical direction are zero flux on the 
vertical boundaries (zero vertical velocity on the plates). One only needs to 
restrict vector fields to be even in y 3 for Vt, V% and odd in y 3 for V 3 , and 
double the cylindrical domain to — h < y 3 < +h. 

We choose V (y) in H S (C), s > 5/2. In [M-N-B-Gj . for the case of "non- 
resonant cylinders", that is, non-resonant a = 2ir/h 1 we have established 
regularity for arbitrarily large finite times for the 3D Euler solutions for O 
large, but finite. Our solutions are not close in any sense to those of the 2D 
or "quasi 2D" Euler and they are characterized by fast oscillations in the e3 
direction, together with a large vortex stretching term 

5V <9V dV 

v(t,y)-W(t,y)=uj 1 —+uj 2 —+oj 3 — , t>0 
oyi oy 2 oy 3 



with leading component 



1. There are no assumptions on 



oscillations in j/i, y 2 for our solutions (nor for the initial condition Vo(j/))- 

Our approach is entirely based on sturying fast singular oscillating limits 
of Eqs. p.l[) - p.5a[) . nonlinear averaging and cancelation of oscillations in the 
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nonlinear interactions for the vorticity field for large fl. This has been devel- 
oped in |B-M-N2j . |B-M-N3j . and [B-M-N4] for the cases of periodic lattice 
domains and the infinite space R 3 . 

It is well known that fully three-dimensional initial conditions with uni- 
formly large vorticity excite fast Poincare vorticity waves |B-M-N2] , [B-M-N3] , 
[B-M-N4 , |Poi] . Since individual Poincare wave modes are related to the 
eigenfunctions of the curl operator, they are exact time-dependent solutions 
of the full nonlinear 3D Euler equations. Of course, their linear superposition 
does not preserve this property. Expanding solutions of (|2.1[) - (|2.8[) along such 
vorticity waves demonstrates potential nonlinear resonances of such waves. 
First recall spectral properties of the curl operator in bounded, connected 
domains: 

Proposition 2.1 ( M-N-B-G ) The curl operator admits a self-adjoint ex- 
tension under the zero flux boundary conditions, with a discrete real spectrum 
X n = ±|A„|, |A„| > for every n and |A n | — * +00 as \n\ — > 00. The corre- 
sponding eigenfunctions 4?^ 



J° = I U e L 2 {D) : V • U = and U • N\ dD = and / U dz = 



Remark 2.2 In cylindrical domains, with cylindrical coordinates (r, 9, z), the 
eigenfunctions admit the representation: 



*&n 1 ,n2,n 3 (^r,ni ,n 2 ,n 3 (f) j *&6,ni ,ri 2 ,n 3 ( r ) j ^z,ni ,n 2 ,n 3 (j") ) 6 6 > (2-15) 



with n2 = 0, ±1, ±2, «3 = ±1, ±2, ... and n\ = 0, 1, 2, .... Here n± indexes 
the eigenvalues of the equivalent Sturm- Liouville problem in the radial coor- 
dinates, and n = (m, 112,113). See |M-N-B-G] for technical details. From now 
on, we use the generic variable z for any vertical (axial) coordinate j/3 or 23. 
For = (vertical averaging along the axis of the cylinder), 2- Dimensional, 
3- component solenoidal fields must be expanded along a complete basis for 
fields derived from 2D stream functions: 



curl*± = ±|A„|*± 



(2.13) 



are complete in the space 





(2.14) 



*„ = ((curl(4) n e 3 ), (j> n e 3 )), <j> n = <f> n (r,8), 
-A(f> n = Hn<t> n , 4> n \dr=0, and 
curl$ n = ({curl((f) n e 3 ), ^„</> n e 3 )). 



Here ((a, 6e 3 )) denotes a 3-component vector whose horizontal projection is 
a and vertical projection is be 3 . 
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Let us explicit elementary swirling wave flows which are exact solutions to 
(ED) and (1271) : 

Lemma 2.3 For every n — (ni,n-2,n^), the following quasiperiodic (T 2 in 
time ) solenoidal fields are exact solution of the full 3D nonlinear Euler equa- 
tions nn\): 



v 0> V) = ^ e 3 x V + exp(^ Ji)*„(exp(-^ Jt)y) exp(±i-^-aClt), (2.16) 

77,3 is the vertical wave number of & n and exp(^Jt) the unitary group of rigid 
body rotations: 





-1 






i 





i) 


, e njt / 2 = I 












cos(^) -sin(f) 
sin(f ) cos(^) 
1 



(2.17) 



Remark 2.4 These fields are exact quasiperiodic, nonaxisymmetric swirling 
flow solutions of the 3D Euler equations. For ^ 0, their second components 

V(t, y) = exp(^ Ji)*„(exp(-^ 3t)y) exp(±^ofit) (2.18) 
are Beltrami flows (curW x V = 0) exact solutions of {2.7]) with V(t = 0, y) — 

V(t,y) in Eq. (|2. 18[) are dispersive waves with frequencies § and p^rfi, 

where a = 3?. Moreover, each V(t, y) is a traveling wave along the cylinder's 
axis, since it contains the factor 

( ^ N 

exp I ian3(±z ± — — -t) 

Note that n.3 large corresponds to small axial (vertical) scales, albeit < 
a|n3/A„| < f. 

Proof of Lemma \2.3\ Through the canonical rigid body transformation for 
both the field V(t, y) and the space coordinates y = (2/1,2/2,2/3): 

V(t, 2/) = e+^Uit, e-^ 2 y) + ^Jy, x = e^^y, (2.19) 
the 3D Euler equations (|2.1| . ()2.2|) transform into: 

d t U + {curlV + fie 3 ) xU = -v|p- ^-(l*i| 2 + M 2 ) + , (2.20) 

V • U = 0, U(t, a;)|t=o = U(0) = V (x), (2.21) 
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For Beltrami flows such that curlXJ x U = 0, these Euler equations (|2.20p ~ 
(I2.2ip in a rotating frame reduce to: 

<9 t U + Oe 3 x U + Vtt = 0, V • U = 0, 

which are identical to the Poincare-Sobolev nonlocal wave equations in the 
cylinder |M-N-B-Oj . [Pol] . [Sob] . |Ar-Khej : 



<9 t * + Qe 3 x * + Vtt = 0, V • * = 0, (2.22) 

d 2 d 2 
W2 curl 2 *-n 2 - 



— curl 2 * - n 2 — * - 0, * ■ N\ dD = 0. (2.23) 



It suffices to verify that the Beltrami flows * n (t, x) = & n (x) exp (^tijy^fltj , 

where $^(x) and ±|A n | are curl eigenfunctions and eigenvalues, are exact 
solutions to the Poincare-Sobolev wave equation, in such a rotating frame of 
reference. 

Remark 2.5 The frequency spectrum of the Poincare vorticity waves (solu- 
tions to (|2.22[) ) is exactly ±ijj^f2, n = (m, U2,n^) indexing the spectrum 
of curl. Note that n 3 = ( zero frequency of rotating waves ) corresponds to 
2- Dimensional, 3-Components solenoidal vector fields. 



We now transform the Cauchy problem for the 3D Euler equations (|2.1[) - 
(|2.2p into an infinite dimensional nonlinear dynamical system by expanding 
V(i, y) along the swirling wave flows (|2.16[) - (|2.18[) : 



V(t,2/) = ~e 3 x y 




+ exp(-Jt)^ ]T u„(t)exp ( ±ijy-^m) *„ ( exp ( -— Jt ) y 

n=(n 1,712,^-3) 



V(< = 0,y) = ^xj + Vofe) 



v„(y) = £ Un(0)*n(y), 

n=(n 1 ,n 2 ,n 3 ) 



(2.24a) 

) 

(2.24b) 
(2.24c) 

(2.24d) 



where denotes the curl eigenfunctions of Proposition 2.1 if 713 ^ 0, and 
*n = ((cwr/(</)„e3), </>„e 3 )) if 713 = (2D case, Remark 2.2). 



As we focus on the case where helicity is conserved for (|2.1[) - (12.2|) . we 
consider the class of initial data Vo such that [M-N-B-G] : 

curlVo • N = on T, 
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where T is the lateral boundary of the cylinder. 

The infinite dimensional dynamical system is then equivalent to the 3D 
Eulcr equations (|2.1[) - (|2.2p in the cylinder, with n = (n\, i%2, TI3) ranging over 
the whole spectrum of curl, e.g.: 



du n ^ /./ n 3 k 3 m 3 

exp 1 ±— — - ± —— ± — — - ailt 



dt p \ \ |A„| A fe |A r 



k3+m.3=n3 
k 2 +m2=n 2 



(2.25) 



x < curl$ k x * m ,*„ > u k (t)u m (t) 



here 



cur/*£ = ±A fe *± if k 3 £ 0, 
curl$ k = ((curl((f> k e 3 ), fi k <p k e 3 )) if k 3 = 

(2D, 3-components, Remark 2.2), similarly for m 3 = and n 3 — 0. The inner 
product < , > denotes the L 2 complex-valued inner product in D. 

This is an infinite dimensional system of coupled equations with quadratic 
nonlinearities, which conserve both the energy 

£(i) = ]T|u ri (i)| 2 

n 

and the helicity 

Hel(t) = Y,±\K\\u±(t)\\ 

n 

The quadratic nonlinearities split into resonant terms where the exponential 
oscillating phase factor in (|2.25p reduces to unity and fast oscillating non- 
resonant terms (f2 >> 1). The resonant set K is defined in terms of vertical 
wavenumbers k 3 ,m 3 ,n 3 and eigenvalues ±\ k , ±A m , ±A„ of curl: 

K = {±-i ± ^ ± ^ = 0, n 3 = k 3 + m 3 , n 2 = k 2 + m 2 }. (2.27) 

*k A m A n 

Here k 2l m2^n 2 are azimuthal wavenumbers. 

We shall call the "resonant Euler equations" the following oo-dimensional 
dynamical system restricted to (k, m, n) € K: 

^ + V < curl$ k x > u fe u m = 0, (2.28a) 

dt 4 

(k,m,n)€iK 

u„(0)=<V ,* n >, (2.28b) 



here curl^ = ±A fc ** if k 3 ^ 0, curl$ k = ((curl((f> k e 3 ), ^ k <p k e 3 )) if k 3 = 0; 
similarly for m 3 — and n 3 = (2D components, Remark 2.2). If there 
are no terms in (|2.28a|) satisfying the resonance conditions, then there will be 
some modes for which — tt = 0. 
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Lemma 2.6 The resonant 3D Euler equations \2.28\) conserve both energy 
E(t) and helicity Hel(t). The energy and helicity are identical to that of the 
full exact 3D Euler equations Ij2.1)) -S2.2 \) . 



The set of resonances K is studied in depth in |M-N-B-G| . To summarize, 
K splits into: 

(i ) 0-wave resonances, with 713 = k 3 = m 3 = 0; the corresponding reso- 
nant equations are identical to the 2-Dimensional, 3-Components Euler 
equations, with I.C. 



1 >" 
h 



Vo(yi,y2,y 3 )dy 3 . 



(ii) Two- Wave resonances, with £37713713 = 0, but two of them are not null; 
the corresponding resonant equations (called "catalytic equations" ) are 
proven to possess an infinite, countable set of new conservation laws 
[M-N-B-Gj . 

(iii) Strict three- wave resonances for a subset K* C K . 
Definition 2.7 The set K* of strict 3 wave resonances is: 

K* = \ ±^- ± ^ ± ^ = 0, k 3 m 3 n 3 ^ 0, n 3 = k 3 + m 3 ,n 2 = k 2 + m 2 

(2.29) 

Note that K* is parameterized by h/R, since a = =£■ parameterizes the eigen- 
values A„, Afc, A TO of the curl operator. 

Proposition 2.8 There exist a countable, non-empty set of parameters for 
which K* 7^ 0. 

Proof. The technical details, together with a more precise statement, are 
postponed to the proof of Lemma 13.71 Concrete examples of resonant ax- 
isymmetric and helical waves are discussed in (MahJ ( cf. Figure 2 in the 
article) . 

Corollary 2.9 Let J Q Vo(yi, y 2 , y 3 ) dy 3 = 0, i.e. zero vertical mean for the 
I.C. V (y) in PO)) . ifP|) , $2.24d\ ) and HIM) . Then the resonant 3D Euler 
equations are invariant on K* : 

~Jf + Yl A fe < * fe x * m , *„ > u fe u m = 0, k 3 m 3 n 3 ^ 0, (2.30a) 

u„(0) =< V ,*„ > (2.30b) 

(where Vp has spectrum restricted to n 3 ^ 0). 



12 



Proof. This is an immediate corollary of the "operator splitting" Theorem 
3.2 in [MZEEQ . ■ 

We shall call the above dynamical systems the "strictly resonant Euler 
system" . This is an oo-dimensional Riccati system which conserves Energy 
and Hclicity. It corresponds to nonlinear interactions depleted on K*. 

How do dynamics of the resonant Euler equations (|2.28[) or (|2.30[) approx- 
imate exact solutions of the Cauchy problem for the full Euler equations in 
strong norms? This is answered by the following theorem, proven in Section 



• Let ~w(t,x) denote the solution to the resonant 3D Euler equations with 
Initial Condition w(0,x) = w(0,j/) = V (y). 

• Let ||w(t,»)|| H - < M S (T M ,M°) on < t < T M , s > 7/2. 
Then, Ve > 0, 3 Q*(T M) M°,e) such that, Vfi > 0*: 



on0<t< T M , V/3 > 1, /3<s-2. Here || • \\ H is defined in FTIS)) . 

The 3D Euler flow preserves the condition curlVo • N = on T, that is 
curlV(t,y) ■ N = on T, for every t > |M-N-B-G| . The proof of this 
"error-shadowing" theorem is delicate, beyond the usual Gronwall differential 
inequalities and involves estimates of oscillating integrals of almost periodic 
functions of time with values in Banach spaces. Its importance lies in that 
solutions of the resonant Euler equations (|2.28[) and/or (|2.30[) are uniformly 
close in strong norms to those of the exact Euler equations (|2.ip - (|2.2j) . on 
any time interval of existence of smooth solutions of the resonant system. 
The infinite dimensional Riccati systems (|2.28[) and (|2.30[) are not just hydro- 
dynamic models, but exact asymptotic limit systems for O ^ 1. This is in 
contrast to all previous literature on conservative 3D hydrodynamic models, 
such as in |G-D-Oj . 



4 of [M-N-B-Gj : 



Theorem 2.10 Consider the initial value problem 
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3 Strictly resonant Euler systems: the SO(3) 
case 



We investigate the structure and the dynamics of the "strictly resonant Euler 
systems" (|2.30|) . Recall that the set of 3- wave resonances is: 

K* = j(fc,m,n) : ± ± ^ = 0, k 3 m 3 n 3 ^ 0, 

Afc A m A n (3-1) 

"3 = k 3 + rn 3 , n 2 = k 2 + to 2 j. 

From the symmetries of the curl eigenfunctions <&„ and eigenvalues A n in the 
cylinder, the following identities hold under the transformation n 2 — ► — n 2 , 

n 3 -> -n 3 

*(ni,-n 2 ,-n 3 ) = **(ni,n 2 ,n 3 ), 
A(m,-n 2 ,-n 3 ) = A(ni,n 2 ,n 3 ) . 

where * designates the complex conjugate (see Section 3, jM-N-B-G] for 
details). The eigenfunctions «&(ni,n 2 ,n 3 ) involve the radial functions 
J n2 (f}(ni,n 2 ,an 3 )r) and J n2 ((3(ni,n 2 ,an 3 )r), with 

A 2 (ni,7i2,n 3 ) = 2 {n u n 2 ,an 3 ) +a 2 n 3 ; 

(3(ni,n 2 ,an 3 ) are discrete, countable roots of equation (3.30) in [M-N-B-G] . 
obtained via an equivalent Sturm-Liouville radial problem. Since the curl 
eigenfunctions are even in r — > — r, ni — > —n\, we will extend the indices 
n\ = 1,2, ...,+oo to — n\ = —1,-2,... with the above radial symmetry in 
mind. 

Corollary 3.1 The 3-wave resonance set K* is invariant under the symme- 
tries o~j, j = 0, 1, 2, 3, where 

cro(ni,n 2 ,n 3 ) = (ni,n 2 ,n 3 ), 
ai(ni,n 2 ,n 3 ) = n 2 , n 3 ), 

a 2 (ni,n 2 ,n 3 ) = (ni,-n 2 ,n 3 ) 
cr 3 (m,n 2 ,n 3 ) = (ni,n 2 ,-n 3 ) . 

Remark 3.2 For < i < 3, < j < 3, < / < 3 a] = Id, a l a ] = -at if 
i 7^ j and o~i<jjO~i — —Id, for i ^ j =/= I. The <jj do 'preserve the convolution 
conditions in K* . 

We choose an a for which the set K* is not empty. We further take the 
hypothesis of a single triple wave resonance (fc, to, n), modulo the symmetries 

(Tf 

Hypothesis 3.3 K* is such that there exists a single triple wave number 
resonance (n,k,m), modulo the symmetries o~j, j — 1,2,3 and crAk) ^ 
k, crj (to) ^ to, (Tj (n) ^ n for j = 2 and j = 3. 
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Under the above hypothesis, one can demonstrate that the strictly resonant 
Euler system splits into three uncoupled systems in C 3 : 

Theorem 3.4 Under hypothesis \3.3l the resonant Euler system reduces to 
three uncoupled rigid body systems in C 3 : 

+ i(Xk — X m )CkmnUkU m = (3.3a) 

dUk 

—± - i(X m - \ n )C krnn U n U; n = (3.3b) 

- %{X n - Xk)Ck mn U n Uk = (3.3c) 

where Ckmn — i < x $m,^Ji >, Gkmn real and the other two uncoupled 
systems obtained with the symmetries o~2{k, m, n) and 03 (fc, to, n). The energy 
and the helicity of each subsystem are conserved: 

j t {U k Ut + U m U* m + U n U*) = Q, 

-ir{XkU k Ut + X m U m U* m + X n U n U*) = 0. 
at 

Proof. It follows from U-k — U£, X(—k) = X(+k), similarly for to and 
n; and in a very essential way from the antisymmetry of < &k x 3> m , >, 
together with curl&k = Xk&k- That Cfc mn is real follows from the eigenfunc- 
tions explicited in Section 3 of [M-N-B-G] . ■ 

Remark 3.5 This deep structure, i.e. SO(3; C) rigid body systems in C 3 is 
a direct consequence of the Lame form of the full 3D Euler equations, cf. Eqs. 
il.3\ ) and \2. 7| ), and the nonlinearity curW x V. 

The system (|3.3[) is equivariant with respect to the symmetry operators 
(21,22,23) -> (2*, 22,23), (2:1,22,23) -» (exp(ixi)zi,exp(ix2)22,exp(ix 3 )2 3 ) , 

provided xi = X2 + X3- It admits other integrals known as the Manley- 
Rowe relations (see, for instance Wc-Wil ). It differs from the usual 3- 
wave resonance systems investigated in the literature, such as in |Zak-ManT] , 
|Zak-Man2j . |Gu-Maj in that 

(1) helicity is conserved, 

(2) dynamics of these resonant systems rigorously "shadow" those of the 
exact 3D Euler equations, see Theorem 12. 101 

Real forms of the system (I3.3[) are found in Gledzer et al. [G-D-Oj . corre- 
sponding to the exact invariant manifold U~k €E iR, U m £ R, U n S R, albeit 
without any rigorous asymptotic justification. The C 3 systems Q3.3P with 
helicity conservation laws are not discussed in [G-D-O] . 
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The only nontrivial Manley-Rowe conservation laws for the resonant sys- 
tem (|3.3[) . rigid body S*0(3; C), which are independent from energy and he- 
licity, are: 

37 (r k r m r n sin(6>„ - 6 k - m )) = 0, 
at 

where Uj = rj exp(i6j), j — k, m, n, and 

£\ = (Aft — A TO )r„ ~ (X m — X n )r k , 
£■2 = (A m - A n )r fc - (A„ - Afc)r m . 



The resonant system (J3T3J) is well known to possess hyperbolic equilibria 
and heteroclinic/homoclinic orbits on the energy surface. We are interested 
in rigorously proving arbitrary large bursts of enstrophy and higher norms 
on arbitrarily small time intervals, for properly chosen h/R. To simplify 
the presentation, we establish the results for the simpler invariant manifold 
U k G «R, and U m , U n £ R. 



Rescale time as: 



Start from the system 



t -> t/C k 



U n + «(A fe - X m )U k U m = 

U k - i(X m - \ n )U n U* m = (3.4) 
U m - i(\n - X k )U n U* k = 

Assume that U k G iR and that U m , U„ G R: set p = iU k , q = U m and r = U n , 
as well as Afe = A, X rn — /i and A„ = v. then 

p + (fj, — v)qr = 

q+{v-X)rp = Q (3.5) 
r + (A — n)pq = 

This system admits two first integrals: 

E = p 2 + q 2 + r 2 (energy) 

9 9 9 , v ( 3 -6) 

H = Xp + fiq 2 + vr 2 (helicity) 

System p.5p is exactly the 50(3, R) rigid body dynamics Euler equations, 
with inertia momenta Ij = j = k,m,n [Arlj . 

Lemma 3.6 ( [Arlj . [G-D-Oj ) With the ordering X k > X m > X n , i.e. X > \i > 
v, the equilibria (0,±1,0) are hyperbolic saddles on the unit energy sphere, 
and the equilibria (±1,0,0), (0,0, ±1) are centers. There exist equivariant 
families of heteroclinic connections between (0, +1,0) and (0,-1,0). Each 
pair of such connections correspond to equivariant homoclinic cycles at (0, 1, 0) 
and (0,-1,0). 
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We investigate bursting dynamics along orbits with large periods, with 
initial conditions close to the hyperbolic point (0, E(0), 0) on the energy sphere 
E. We choose resonant triads such that Afe > 0, A„ < 0, \ k ~ |A n |, |A m | -C X k , 
cquivalently: 

A > // > v, \v < 0, \n\ < A and A ~ \v\. (3.7) 

Lemma 3.7 There exist h/R with K* ^ 0, such that 

Afc > A m > A„, A fe A„ < 0, |A m | < A fc and A fc ~ |A n |. 

Remark 3.8 Together with the polarity ± of the curl eigenvalues, these are 
3-wave resonances where two of the eigenvalues are much larger in mod- 
uli than the third one. In the limit \k\, \m\, \n\ 3> l,Afe ~ ± j fc| , A m ~ 
±|m|, A„ ~ ±M, the eigen] 'unctions 3? have leading asymptotic terms which 
involve cosines and sines periodic in r, cf. Section 3 [M-N-B-G] . In the 
strictly resonant equations \2.30\) . the summation over the quadratic terms 
becomes an asymptotic convolution in n\ = k\ +n\. The resonant three waves 
in Lemma \3. 7| are equivalent to Fourier triads k + m = n, with \k\ ~ \n\ 
and \m\ -C \k\, \n\, in periodic lattices. In the physics of spectral theory of 
turbulence [Frij , [Lesj , these are exactly the triads responsible from transfer of 
energy between large scales and small scales. These are the triads which have 
hampered mathematical efforts at proving the global regularity of the Cauchy 
problem for 3D Navier-Stokes equations in periodic lattices Fe . 

Proof of Lemma\3~T\ QM-N-B-Gj ) The transcendental disp ersion law for 
3- waves in K* for cylindrical domains, is a polynomial of degree four in #3 = 
l/h 2 : 

P(tf 3 ) = P«& + P 3 4 + P2-&1 + Pl#3 + Po = 0, (3.8) 
with n-i = &2 + rri2 and 713 = k% + 7713. 

Then With h k = tlhl^hsl hm = P 2 (mi,m2, a m 3 ) h = /3 2 (ni,» 3 ,an 3 ) f 

the radial Sturm-Liouville problem in Section 3, [M-N-B-G] . the coefficients 
of P("&3) are given by: 

Pi = -3, 

P3 = -4(hk + h m + h n ), 

P 2 = -6(hkh m + hkh n + h m h n ), 

Pi — —12hkh m h n , 

Po = h 2 m h 2 n + h 2 k h 2 n + h 2 m h\ - 2{h k h m h 2 n + h k h n h 2 m + h rn h n h 2 k ). 

Similar formulas for the periodic lattice domain were first derived in [B-M-N2] , 
[B-M-N3 , B-M-N4 . In cylindrical domains the resonance condition for K* 
is identical to 



05>3 + h k v% + h m VW+h 
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with tf 3 = h k = P 2 {k)/kl h m = 2 {m)/m 2 3 , h n = (3 2 {n)/n 2 - Eq. (Ml is 
the equivalent rational form. 

From the asymptotic formula (3.44) in [M-N-B-G] . for large (3: 

IT 7T 

/3(n 1 ,n 2 ,n 3 ) ~ ni7r + n 2 - + - + V) (3-9) 

where tp = if lim^^ = (e.g. h fixed, m 2 /m 3 — > 0) and i/j = ±-| if 

lim^-^r = ±oo (e.g. ^| fixed, /i — > oo). The proof is completed by taking 

leading terms P + $ 3 Pi in |3T8|) . # 3 = < 1, and m 2 = 0, fc 2 = 0(1), n 2 = 
0(1). ■ 

We now state a theorem for bursting of the H 3 norm in arbitrarily small 
times, for initial data close to the hyperbolic point (0, E(0), 0): 

Theorem 3.9 (Bursting dynamics in H 3 ). Let A > [i > v, Xv < 0, <C A 
and A - Let W(t) = X 6 p(t) 2 + ^q{t) 2 + v & r{t) 2 the H 3 -norm squared of 
an orbit of L3.5\) . Choose initial data such that: W(0) = A 6 p(0) 2 + /i 6 g(0) 2 
with X 6 p(0) 2 ~ \W{Q) and fi 6 q(0) 2 ~ \W{Q). Then there exists t* > 0, such 
that 

W (t)>UjA w(o) 

where t* < Ln{X/\^\){XM)- 1 . 

Remark 3.10 Under the conditions of Lemma \3.1I\ {j^j ^ 1> whereas 

/i 2 (Ln(A/|^|))(A/|/i|) _1 -C 1. Therefore, over a small time interval of length 
0(^ 2 (Ln(A/|^|))(A/|/x|)- 1 ) < 1, ifte ratio ||U(t)|| H 3/||U(0)|| ff 3 grows up to 
a maximal value O ((A/|/i|) 3 ) S> 1. Since the orbit is periodic, the H 3 semi- 
norm eventually relaxes to its initial state after some time (this being a mani- 
festation of the time-reversibility of the Euler flow on the energy sphere). The 
"shadowing" theorem \2.10\ with s > 7/2 ensures that the full, original 3D Eu- 
ler dynamics, with the same initial conditions, will undergo the same type of 
burst. Notice that, with the definition il.13]) of \\ ■ \\h b > °ne has 

I |Oe 3 x y\\ H a = \\curl 3 (fle 3 x y)\\ L i = 0. 

Hence the solid rotation part of the original 3D Euler solution does not con- 
tribute to the ratio \ \V(t)\ \ H3 /\ |V(0)| \ H s . 



Theorem 3.11 (Bursting dynamics of the enstrophy). Under the same con- 
ditions for the 3-wave resonance, let 3(t) = X 2 p(t) 2 + fi 2 q(t) 2 + v 2 r(t) 2 the 
enstrophy. Choose initial data such that 3(0) = A 2 p(0) 2 + /i 2 g(0) 2 + i/ 2 r(0) 2 
with X 2 p(0) 2 ~ §3(0), Ai 2 g(0) 2 - §3(0). Then there exists t** > 0, such that 

'A 



3 (n >. /( 



where t** < ^-7=Ln (A/| M |) (A/|/x|) 



-l 
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Remark 3.12 It is interesting to compare this mechanism for bursts with ear- 
lier results in the same direction obtained by DiPerna and Lions. Indeed, for 
each p € (1, oo), each 5 € (0, 1) and each t > 0, Di Perna and Lions [DiPe-Lij 
constructed examples of 2D-3 components solutions to Euler equations such 
that 

l|V(0)||v»ri* < e while \\V(t)\\ w i, P > 1/5. 
Their examples essentially correspond to shear flows of the form 

I u{x 2 ) 
V{t, Xl ,x 2 ) = 

\w(xi - tu(x 2 ),x 2 ) / 

where u £ W^f while w € W^f . Obviously 

((d 2 — tu'(x 2 )di)w(xi — tu{x 2 ),x 2 ) 
—diw(xi — tu(x 2 ), x 2 ) 
-u'{x 2 ) 

Thus, all components in curl~V(t,Xi,x 2 ) belong to L p oc , except for the term 

-tu 1 \x 2 )d\w{xi - tu{x 2 ), x 2 ) . 

For each t > 0, this term belongs to LP for all choices of the functions u G 
W^ p and w £ Wf.f X2 if and only if p — oo. Whenever p < oo, DiPerna and 
Lions construct their examples as some smooth approximation of the situation 
above in the strong W 1,p topology. 

In other words, the DiPerna-Lions construction works only in cases where 
the initial vorticity does not belong to an algebra — specifically to LP , which 
is not an algebra unless p — oo. 

The type of burst obtained in our construction above is different: in that 
case, the original vorticity belongs to the Sobolev space H 2 , which is an algebra 
in space dimension 3. Similar phenomena are observed in all Sobolev spaces 
with (3 > 2 — which are also algebras in space dimension 3. 

In other words, our results complement those of DiPerna-Lions on bursts 
in higher order Sobolev spaces, however at the expense of using more intricate 
dynamics. 

We proceed to the proofs of Theorem 13.91 and 13.111 We are interested in 
the evolution of 

5 = X 2 p 2 + [i 2 q 2 + v 2 r 2 (enstrophy) (3.10) 

Compute 

S = -2 (A 2 (Ai - v) + ii 2 (y - A) + v 2 (\ - fi)) pgr (3.11) 

then 

(pgr) = -(p - v)g 2 r 2 - (y - X)r 2 p 2 - (A - fi)p 2 q 2 (3.12) 
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Using the first integrals above, one has 

Van q 2 
\r 2 

where Van is the Vandermonde matrix 

Van - 



(3.13) 



1 1 1 

A fx v 
A 2 u 2 v 2 



For A 7^ /U 7^ v ^ A, this matrix is invertible and 
Van" 1 = 















1 




(A- 


M)(A- 
i/A 




(A- 


u)(A-i/) 
("+A) 


(A- 


-m)(a- 
1 








-A) 




")(M-A) 
(A+m) 


(M- 




A) 




A)(„- 


/*) 




A)(i/-^) 




A)(„- 


/<) 



Hence 



P = 



so that 



\ 
/ 

(S - (li + v)H + llvE) 
{E-{v + X)H + vXE) 

(H - (i/ + A)if + i/Af?) (H - (A + n)H + XllE) 



(A- 


/*)(A-i/) 




1 


(/x- 


!/)(/* -A) 




1 



(3.14) 



(A - 


li)(X - v)(li - 






(S-(A + M )iJ + 


- Aui?) (S - (u 


+ v)H + 


LlvE) 


(A- 


li){X-v)(li- 


v) 




(H-( M + i/)H-F 


- iivE) (3 - (i/ 


+ X)H + 


v\E) 


(A - 


m)(a - /^)(m - 







(li - z/)q 2 r 2 
(i/ - A)r 2 p 2 
(A - u)pV 

Later on, we shall use the notations 

£- (A, u, v) = (fi + i>)H — livE 

x (A, a, v) = (u + X)H - liXE (3.15) 
x + (\, a, v) = (A + v)H - XvE 

Therefore, we find that 5 satisfies the second order ODE 

5 = - 2Kx tfl , v ((H - X-(X,li, u))(E - x (X,li, v)) 

+(E - x (A, /x, i/))(S - .t+(A, u, i/)) + (S - x + (X, li, z/))(S - z (A, u, i/))) 

which can be put in the form 

S = -2K X ^ V P' X ^ V {~) (3.16) 
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where P\^m is the cubic 

Px.^AX) = (X- x-(\, fi, u))(X - x {\, /j,, v)){X - x+(X, ft, u)) (3.17) 
and 

X 2 {p-v)+p 2 {v-X) + v 2 {X- l ,) 

K \,li,v = 7T T7T w x (3.18) 

In the sequel, we assume that the initial data for (p, q, r) is such that 
r(0)=0, P(0)(g(0)^0 

Let us compute 

x- (A, n, v) = Xup(0) 2 + fi 2 q{0) 2 + fj,(X - v)p{0) 2 
x (A, fi, v) = A 2 p(0) 2 + p 2 q(0) 2 



x + (X, M , «/) = A 2 P (0) 2 + ( ^ - ^ ) A(0) 2 

A 4 M 



(3.19) 



We shall also assume that 

X>n>v, Xv<0, \fj,\ <C A and A ~ \v\ (3.20) 

Then K\ tlJ , t v > — in fact K\ jfl u ~ 2, and 5 is a periodic function of i such 
that 

inf E(t) = xo(X, fi, v) , supS(i) = x + (X, fi, v) (3.21) 
with half-period 

Ta -- ^ wtT= / / p ^ (3 ' 22) 

We are interested in the growth of the (squared) H 3 norm 

W{t) = X 6 p(t) 2 + n\{t) 2 + is 6 r(t) 2 (3.23) 

Expressing p 2 , q 2 and r 2 in terms of E, H and 5, it is found that 

= X 6 (E.-x-(X,fj.,i/)) M 6 (S~a: + (A,/x,;/)) ;/ 6 (S - gp(A, /mQ) 
' (A-/XKA-!/) + (/x-i/JOi-A) + {v-X^y-v) l ' J 

Hence, when S = x + (X, fi, v), then 

= A 6 (x+(A,/x^) -x_(A,/x,t/)) ;/ 6 (x + (A,^^) - x (A,//, ^)) 
(A-mKA-i/) + [y — X){y — /j,) 

X Q {x+{X,nv) -x_{X,n,v)) 
(A-/x)(A-i/) 
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> 



Let us compute 

x + (A, \iv) — X- (A, v,v) = (A - /x)(A - ^(0) 2 + f — - — - 1 ) A(0) 2 
> -^(0) 2 ~ A 2 q(0) 2 

(3.25) 

We shall pick the initial data such that 

W(0) = A 6 p(0) 6 + A(0) 6 with A 6 p(0) 2 ~ |W(0) and A(0) 2 ~ |W(0) 

(3.26) 

Hence, when S reaches a;+(A, /u, i/), one has 

W S (A-ri(A-,) ~ ^(X-MX-,) Wm ~ ^ H ' (0) ■ <3 ' 27) 

Hence jumps from W(0) to a quantity ~ jj^W(O) in an interval of time 
that does not exceed one period of the S motion, i.e. 2T\ tflt „. Let us estimate 
this interval of time. We recall the asymptotic equivalent for the period of an 
elliptic integral in the modulus 1 limit. 

Lemma 3.13 Assume that x_ < x n < x+. Then 
dx 1 



In 



Eo y/(x-X-)(x-X )(x+-x) V X + - X - 

uniformly in X-, xq, and x + as — ► 1. 
Here 



X + — Xq 
X + —X- 



i < i ImI 3 



Next 
so that 



y/x + (\,H,v)-X-(\,tJL,v) ~ V^MOF A V W ^ 

x (A, n, v) - x-(X, n, v) = (X- ;t) (A - v)p(0) 2 (3.28) 



, _ x+-xo , _ /, (\-n)(\-v)p(0) 2 

" x+-x- 1 V 1 (\-^)(x-u)p(ay i +{^(u+X)- v \-^)q{^ 

(A - /i)(A - u)p(0) 2 + {j4y + A) - i/A - M 2 )g(0) 

2(A- M )(A-^(0) 2 
g(0) 2 iW^OVV _ 
~ 2p(0) 2 ~ 2 W(0)/2A 6 ~ V 
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Hence 

^^wwM^^ww,^ 1 ^ (3 - 29) 



Conclusion: collecting (|3.26|) , (|3.27p and (|3.29[) . we see that the squared H 3 
norm W varies from W(0) to a quantity ~ p e W{0) in an interval of time 

< 12 M 2 lnp^ (Here = A / N 

We now proceed to obtain similar bursting estimates for the enstrophy. 
We return to (|3.21[) and (|3.22[) . Pick the initial data so that 

3(0) = X 2 p(0) 2 + pi 2 q{0) 2 with X 2 p{0) 2 ~ -3(0) and p 2 q(0) 2 ~ ^5(0). 
Then 

x+(X,fi,u) - x-(X,fi,u) 

= (A - (x)(X v) P (0) 2 + (^±* - ^ - l) p 2 q(0) 2 

- 2A 2 p(0) 2 + A 2 c?(0) 2 ~ \s(0) 
M 

while 

x (A, M , i>) - x_ (A, M , v) = (A - ,j) (A - ^(O) 2 ~ 2A 2 p(0) 2 ~ 3(0). 
Hence, in the limit as p = X/\p\ — ► +oo, one has 
2T 11, 1 



2V2 v /^3(o) i-yrr^L 

11 1 1 lnp 

m - 



2723(0)/9 1 - v/1 - 2p- 2 v/23(0) P 
And 3 varies from 

xq(X,h,v) = 3(0) to a:+(A, p, v) ~ p 2 5(0) 
on an interval of time of length T\^ v . I 

4 Strictly resonant Euler systems: the case of 
3-waves resonances on small-scales 

4.1 Infinite dimensional uncoupled SO (3) systems 

In this section, we consider the 3-wave resonant set K* when 

1 

if 



|fc| 2 , H 2 , |n| 2 > — , <? 7«1, 
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i.e. 3-wave resonances on small scales; here \k\ 2 = k\ + fc| + fcf, where 
(k\,k2,k 3 ) index the curl eigenvalues, and similarly for \m\ 2 , \n\ 2 . Recall 
that k2 + m,2 — 77,2, k% + 7773 = 773 (exact convolutions), but that the sum- 
mation on hi, 7711 on the right hand side of Eqs. (|2.30[) is not a convolution. 
However, for \k\ 2 , \m\ 2 , \n\ 2 > -j, the summation in k\, m\ becomes an 
asymptotic convolution. First: 

Proposition 4.1 The set K* restricted to |fc| 2 , |m| 2 , |t7| 2 > V77, < 77 <C 1 
is not empty: there exist at least one h/R with resonant three waves satisfying 
the above small scales condition. 

Proof. We follow the algebra of the exact transcendental dispersion law 
(jni§ derived in the proof of Lemma |3~T1 Note that P($ 3 ) < for tf 3 = ^ 
large enough. We can choose h m = 13 i m ^'^,am 3 ) _ g& y ^ ^ e specific limit 
2 ^ — ► 0, and /3(mi, m 2 , am 3 ) ~ 777 i7r + 7772f + f . Then P a = h? k h 2 n > and 
P(i?3) must possess at least one (transcendental) root $ 3 = ■ 

In the above context, the radial components of the curl eigenfunctions in- 
volve cosines and sines in ^ (cf. Section 3, [M-N-B-G] ) and the summation in 
k\ , mi on the right hand side of the resonant Euler equations (|2.30p becomes 
an asymptotic convolution. The rigorous asymptotic convolution estimates are 
highly technical and detailed in |Fro-M-N] . The 3-wave resonant systems for 
\k\ 2 , |m| 2 , |n| 2 > ^ are equivalent to those of an equivalent periodic lattice 
[0, 2tt] x [0, 2tt] x [0, 2irh], $3 = the resonant three wave relation becomes: 

±(03+01^ +02^ '=0, (4.1a) 

V m 3 m J 

k + 771 = 77, k^m^n^ 7^ 0. (4-lb) 

The algebraic geometry of these rational 3-wave resonance equations has been 
investigated in depth in [B-M-N3j and [B-M-N4] . Here #1, $2, $3 are periodic 
lattice parameters; in the small-scales cylindrical case, x9i = ^2 = 1 (after 
rescaling of n 2 , k 2 , 777 2 ), $3 = l/^ 2 , h height. Based on the algebraic 
geometry of "resonance curves" in [B-M-N3 , |B-M-N4] . we investigate the 
resonant 3D Euler equations (|2.30[) in the equivalent periodic lattices. 

First, triplets (fc, m, n) solution of (|4.ip are invariant under the reflec- 
tion symmetries <7o , a\ , 02 , 03 defined in Corollary 13.11 and Remark 13.21 <jq = 
Id, o-j(k) = (eijki), 1 < i < 3, e^- = +1 if i ^ j, e it j = -1 if i = j, 1 < j < 3. 
Second the set K* in (|4.1[) is invariant under the homothetic transformations: 

(fc, 777, 77) — > (7A;, 7m, 777), 7 rational. (4-2) 
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The resonant triplets lie on projective lines in the wavenumber space, with 
equivariance under cr,- , < j < 3 and 7-rescaling. For every given equivariant 
family of such projective lines, the resonant curve is the graph of ^ versus 
^2-, for parametric domain resonances in $1, $2, $3- 

Lemma 4.2 (p. 17, [B-M-N4] ). For every equivariant (k,m,ri), i/ie resonant 
curve in the quadrant dx > 0, $2 > 0, $3 > is i/ie graph of a smooth function 
$3/$! = F{ , &2/^i) intersected with the quadrant. 

Theorem 4.3 (p. 19, [B-M-N4J. A resonant curve in the quadrant #3/1?! > 0, 
•&2/'&i > is called irreducible if: 





k'2 








:0 














2 


»3 







dei mg TO2 m i ^ °- ( 4 - 3 ) 



irreducible resonant curve is uniquely characterized by six non-negative 
algebraic invariants V\, Vi, 72-i, 7?-2> <Si, 1S2; smc/i £/ia£ 

{$f} = {^f}, {||} = W.^}, 

and permutations thereof. 

Lemma 4.4 (p. 25, |B-M-N4] ). For resonant triplets (k, m, n) associated to 
a given irreducible resonant curve, that is verifying Eq. |^73p ; consider the 
convolution equation n = k + m. Let o"j(n) 7^ n, Vi, 1 < i < 3. TTien f/iere 
are no more that two solutions (k,m) and (m,k), for a given n, provided 
the six non- degeneracy conditions (3.39)-(3.44) in [B-M-N4] for the algebraic 
invariants of the irreducible curve are verified. 

For more details on the technical non-degeneracy conditions, see the Ap- 
pendix. An exhaustive algebraic geometric investigation of all solutions to 
n = k + m on irreducible resonant curves is found in [B-M-N4] . The essence 
of the above lemma lies in that given such an irreducible, "non-degenerate" 
triplet (k, m, n) on K*, all other triplets on the same irreducible resonant 
curves are exhaustively given by the equivariant projective lines: 

(k,m,n) — * (7/0,7771,771), for some 7 rational , (4-4) 
(k,m,n) — > (ajk,crjm,o-jn), j = 1,2,3, (4-5) 

and permutations of k and m in the above. Of course the homothety 7 and 
the o~j symmetries preserve the convolution. This context of irreducible, "non- 
degenerate" resonant curves yields an infinite dimensional, uncoupled system 
of rigid body .SO^jR) and SO(3; C) dynamics for the 3D resonant Euler 
equations (|2.30[) . 
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Theorem 4.5 For any irreducible triplet (k, m, n) which satisfy Theorem \4.3\ 
and under the "non-degeneracy" conditions of ' Lemma \4~4\ (cf. Appendix), the 
resonant Euler equations split into the infinite, countable sequence of uncou- 
pled 50(3; R) systems: 

flfc = rfc m „(A m — A„)a m a„, (4.6a) 
a m = Tkmn(K - Xk)a n ak, (4.6b) 
a n = rfe mn (Afe — A m )afcd m , (4.6c) 

for all (k, m, n) = r y(aj(k*), cr, (m*), aj (n*)), 7 = ±1, ±2, ±3..., < j < 3. 

k* ,m* ,n* are some relatively prime integer vectors in Z 3 characterizing the 
equivariant family of projective lines (k,m,n); Tk mn = i < x <& m ,<I>* >, 
Tfcmn real. 

Proof. Theorem 14.51 is a simpler version for invariant manifolds of more 
general 50(3; C) systems. It is a straightforward corollary of Proposition 3.2, 
Proposition 3.3, Theorem 3.3, Theorem 3.4 and Theorem 3.5 in [B-M-N4] . 
The latter article did not explicit the resonant equations and did not use the 
curl-helicity algebra fundamentally underlying this present work. Rigorously 
asymptotic infinite countable sequences of uncoupled 50(3; R), 50(3; C) sys- 
tems are not derived via the usual harmonic analysis tools of Fourier modes, 
in the 3D Euler context. Polarization of curl eigenvalues and eigenfunctions 
and helicity play an essential role. 

Corollary 4.6 Under the conditions A n * — Afc* > 0, Afc* — A m * > 0, the 
resonant Euler systems J^.6p admit a disjoint, countable family of homoclinic 
cycles. Moreover, under the conditions A„* 3> +1, A m * -C — 1, |Afc*| <C A n », 
each subsystem d4-6| ) possesses orbits whose H s norms, s > 1, burst arbitrarily 
large in arbitrarily small times. 

Remark 4.7 One can prove that there exists some T max , < T rnax < 00, 
such that |rfe mn | < T max , for all (fc,m,n) on the equivariant projective lines 
defined by J^. 7| ). Systems |^.g[ ) "freeze" cascades of energy; their total enstro- 
phy E(t) = E(fc,m,n)( A fc a fcW + A m<4M + remains bounded, albeit 

with large bursts o/S(t)/H(0), on the reversible orbits topologically close to 
the homoclinic cycles. 

4.2 Coupled SO (3) rigid body resonant systems 

We now derive a new resonant Euler system which couples two 50(3; R) 
rigid bodies via a common principle axis of inertia and a common moment 
of inertia. This 5-dimensional system conserves energy, helicity, and is rather 
interesting in that dynamics on its homoclinic manifolds show bursting cas- 
cades of enstrophy to the smallest scale in the resonant set. We consider the 
equivalent periodic lattice geometry under the conditions of Proposition ^. II 
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In Appendix, we prove that for an "irreducible" 3-wave resonant set which 
now satisfies the algebraic "degeneracy" (|A-4[) , there exist exactly two "prim- 
itive" resonant triplets (k,m,n) and (k,m,n), where k, to, k, in are relative 
prime integer valued vectors in Z 3 : 

Lemma 4.8 Under the algebraic degeneracy condition {A-$ the irreducible 
equivariant family of projective lines in K* is exactly generated by the follow- 
ing two "primitive" triplets: 

n = k + to, k = ak, to = bm, (4.8a) 

n = k + m, k = a <7j( k) + b o~j(rn), (4.8b) 

that is, 

n = ak + bm, (4.8c) 

n = a'crj(fc) + b'o-j(m), (4.8d) 

where o~i =^ o~j are some reflection symmetries, a, b, a', b' are relatively 
prime integers, positive or negative, and k, to are relatively prime integer 
valued vectors in Z 3 , that is: 

(a, a') = (6, b') = (a, b) = (a', b') = 1, (k, to) = 1, 

where ( , ) denotes the Greatest Common Denominator of two integers. All 
other resonant wave number triplets are generated by the group actions o~i, 
I = 1,2,3, and homothetic rescalings (k,m,n) — > 7(fc,m,n), [k,rh, n) — » 
7(fc, rh, n), (7 € Z) of the "primitive" triplets. 

Remark 4.9 It can be proven that the set of such coupled "primitive" triplets 
is not empty on the periodic lattice. The algebraic irreducibility condition of 
Lemma \4.2\ implies that ±k3/\k\ =±&3/|fc| and ±TO3/|to| = ±7713/ |m|, which 
is obviously verified in equations f^.#| ). 

Theorem 4.10 Under conditions of Lemma \4-S\ the resonant Euler system 
reduces to a system of two rigid bodies coupled via a n (t): 

(4.9a) 
(4.9b) 

AriOfapm (4.9c) 
(4.9d) 
(4.9e) 

where T — i < <frfc x <& m , 3>* >, f = i < x <fr TO , <&* >. Energy and Helicity 
are conserved. 







An)r^m^n 


dm 


= (A„- 


- Xk)Ta n ak 


a n 


= (A fe - 


- A m )rafca m + (Aj 


Qi-fti 


- (A„ - 


- \%)fa n a~ k 
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= (Am 
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Theorem 4.11 The resonant system possesses three independent con- 

servation laws: 

£ l =a 2 k + {l-a)a 2 m , (4.10a) 

£2 = 4 + « a m + (1 - «)<4. (4.10b) 

£ 3 = a? + ac4, (4.10c) 

where 

a = (A m - A fc )/(A„ - A fe ), (4.11a) 

5 = (A* - A„)/(A S - An). (4.11b) 

Theorem 4.12 Under the conditions 

A m < A fc < A„, (4.12a) 
Am < A„ < Aj., (4.12b) 

which imply a < 0, a < 0, £/ie equilibria (±a.fc(0), 0, 0, 0, ±a^,(0)) are hyper- 
bolic for |«r (0) I small enough with respect to |afc(0)|. The unstable manifolds 
of these equilibria are one dimensional, and the nonlinear dynamics of system 
|^.ff[ ) are constrained on the ellipse £\ f^. 10a\ ) for ak(t), a m (t), the hyperbola 
£3 H4-10c\ ) for a%(t), a m (t), and the hyperboloid £2 (4-10b\ ) for a m {t), a m (t), 
a n (t). 

Theorem 4.13 Let the 2-manifold £\ n £2 H £3 fee coordinatized by (a m , a m )- 
On this 2-manifold, the resonant system is Hamiltonian, and therefore 

integrable. Its Hamiltonian vector field h is defined by 

L h co = r(An-A fe )^-f(A n -A £ )^, (4.13) 

where l^lo designates the inner product of the symplectic 2-form 

dam A da m 



aka ri a^ 

with the vector field h. 



(4.14) 



Proof of Theorem S\ Eliminating dfe(t) via £\, a n (t) via £2, a^(t) via 
£3, the resonant system f)4.9[) reduces to: 

a m = ±r(A n - Afe)(£i - (1 - a)o4j'(£ 2 - ota 2 m + (d - 1)<4J* 
dm = ±f(\ n -\k)(£ 2 -aa m + (a-l)a m )2(£ 3 -aaf h )2-, 

after changing the time variable into 



(£1 - (1 - a)a m ) 2 (£ 2 - aa 2 m + (d - l)a|J « (£3 - fi<4) 2 ds . 



/o 
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On each component of the manifold £ \ PI £ % PI £ 3 , the following functional 
are conserved: 



Observe that the system of two coupled rigid bodies (|4. 9[) does not seem to ad- 
mit a simple Lie-Poisson bracket in the original variables a m , a n , a m , a^). 
Yet, when restricted to the 2-manifold £\ r\£i that is invariant under the 
flow of (|4.9p . it is Hamiltonian and therefore integrable. 

This raises the following interesting issue: according to the shadowing 
Theorem 12.101 the Euler dynamics remains asymptotically close to that of 
chains of coupled 50(3; R) and 5*0(3; C) rigid body systems. Perhaps some 
new information could be obtained in this way. We are currently investigating 
this question and will report on it in a forthcoming publication |G-M-Nj . 

Already the simple 5-dimensional system (|4.9j) has interesting dynamical 
properties, wich we could not find in the existing literature on systems related 
to spinning tops. 

Consider for instance the dynamics of the resonant system (|4.9p with I.C. 
topologically close to the hyperbola equilibria (±0^(0), 0, 0, 0, ±a^(0)). Un- 
der the conditions of (|4.12p and with the help of the integrability Theorem 
14.131 it is easy to construct equivariant families of homoclinic cycles at these 
hyperbolic critical points: 

Corollary 4.14 The hyperbolic critical points (±a/c(0), 0, 0, 0, ±a^.(0)) pos- 
sess 1- dimensional homoclinic cycles on the cones 



with a < 0, a < 0. 

Note that these are genuine homoclinic cycles, NOT sums of heteroclinic 
connections. Initial conditions for the resonant system (|4.9p are now chosen 
in a small neighborhood of these hyperbolic critical points, the corresponding 
orbits are topologically close to these cycles. With the ordering: 




a n + (1 - <S)<4 



2 



(4.15) 




(4.16a) 
(4.16b) 
(4.16c) 
(4.16d) 
(4.16e) 
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which can be realized with | — | ^> 1 and \k-\ <C 1 in the resonant triplets 



(|4.8[) . we can demonstrate bursting dynamics akin to Thcorcm l3.9l and l3.11l for 
cnstrophy and H s norms, s > 2. The interesting feature is the maximization 
of |a£(f)| near the turning points of the homoclinic cycles on the cones (|4.15|> . 
This corresponds to transfer of energy to the smallest scale k, X^. 

In a publication in preparation, we investigate infinite systems of the cou- 
pled rigid bodies equations (|4.9[) . 



APPENDIX 

We focus on a resonant wave number triplet (n, k, m) € (Z*) 3 verifying 
• the convolution relation 

n = k + m, (A-l) 



the resonant 3-wave resonance relation 



n 3 

+ J + 



1,3 v -r u 2 ^ 2 -r v 6 n, 3 

m 3 

± = 0, 

^J-d\m\ + ■d2m\ + $37713 



(A-2) 



'3 

• the condition of "non-catalyticity" 

k 3 m 3 n 3 ^ 0, (A-3) 

• and the degeneracy condition of [B-M-N4] (see p26) 

G l {j{k,m) = kiUjini + kirrijUi — 0, (A-4) 
where I) is a permutation of (1, 2, 3). 

Then, we know (see lemma 3.5 (2) of [B-M-N4] ) that the system of equations 
(|A-3|) - (IA-4|) for the unknown k and m, given the vector n, admits exactly 4 
solutions in Z 3 x Z 3 : 

(k,m), (m,k), (k,m), (rh,k). 
Here k and m are the two vectors of the original resonant triplet, whereas 

k = a<Ji{k), rh — f3aj(m) 

where 

miki — mtki , r , , „ mikj — m,-fc; , r 
a = V_u T ^ * ' ±X > and ^ = ITT V t {°« ±X } 

TOifc; + m;fci TTlikj + ITljki 
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and where the symmetries Ui and Uj are defined by 

a, : u = (M()i=l,2,3 ((— l)* w Ui) i=lj2j3 ' 
One verifies that 

CTj = cr| = Id, <7j(7j = (Jj(Ji = —a\. 

That is, the group generated by Uj and <Tj is the Klein group Z/2Z x Z/2Z. 

Let us first write the irrational numbers a and /3 under the irreducible 
representation 

a' b' 

a = — , /3 = — , with a, a', 6, 6' S Z* and (a, a') = (6, b') = 1, 
a o 

where ( , ) denotes the Greatest Common Denominator of the integer pair. 
From k S Z 3 , it follows that a\a'k; but since (a, a') = 1, the Euclid's lemma 
yields that a\k. Similarly, b\m. Now set 

-In 1 o 

k = -k € Z , m = -m e Z J . 
a o 

Hence the integer vector n admits the two decompositions 
n = ak + bm = affj(fc) + 6'oj (m). 

Since the function 

Z3 



is homogeneous of degree 0, we see that within the resonance condition (|A-2|) 
we can replace each vector k, m and n by any colinear vectors - either integer 
or not. Suppose now that there exists some positive integer d ^ 1 such that 
d\k; then d\n, so that by setting 

n = -n, k = -k, m = -m 
d d d 

we finally obtain 

no = ako + bko = a'a^ko) + b'aj(nio). 

The triplets (no, ako, bmo) and (no, a'ai(ko), b'aj(mo)) further verify from the 
above remark, the convolution relation (|A-1|) and the resonance relation (|A-2j) . 
Hence, without loss of generality, we can assume that the only positive integer 
d such that d\k and d\fn is 1; which we denote by 

(k,rn) = 1. 

Equivalently, 

fciZ + k 2 Z + k 3 Z + m 1 Z + m 2 Z + m 3 Z+ = Z. 
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Finally, suppose there exists some positive integer d ^ 1 such that d\a and 
d\b. Then d\n; set 

n = -n, a Q = -a, b Q = -b. 
d d d 

Observe that 

G-j (a fc, b m) = —Glj(ak,bm) = 0. 

It follows from lemma 3.5 (2) of [B-M-N4] that the vector no of the resonant 
triplet (no, a$k, born) can also be written as 

no = k + m with (no, k, m) verifying (|A-2|) . 

But then 

n = driQ = ak + bm = a <Xj(fc) + b'aj (m) = dk + dm. 

From lemma 3.5 (2) of [B-M-N4] . (dk,drh) must coincide with either one of 
the pairs 

{a! di(k), b' <jj(jn)), (b'(Tj(m), a <Zj(A;)). 

In particular, d\a'k and d\b'rn. Since d\a and (a, a') = 1, we have (d, a'); 
similarly (d, b') = 1. But then Euclid's lemma yields that d\k and d|m, which 
contradicts the fact that (fc,m) = 1. Hence we have proven that (a, 6) = 1. 
In a similar way, one can show that (a', b') = 1. 

Conclusion: It follows from the above study that n G Z* admits the two 
decompositions 

n = ak + bm = a <Ji(k) + b'aj(rn) 

with 

(a, a') = (6, 6') = (a, 6) = (a', b') = 1, (fc, m) = 1. 

The triplets (n,ak,bm) and (n,aai(k),b'aj(m)) both verify the resonant 
condition (|A-2|) (from the homogeneity of this condition) as well as the condi- 
tion of non-catalyticity (|A-3[1 . Indeed, aba'b' ^ and the condition (|A-3|) on 
the initial triplet (n, k, m) imply that the reduced triplet (n, k, m) also verifies 
(|A-3jl ). Finally, the degeneracy condition (|A-4|) 

Gf tj {ak,bm) = 

is verified. 
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